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Abstract 

The classical flexure problem of non- linear incompressible elasticity 
is revisited for elastic materials whose mechanical response is different 
in tension and compression — the so-called bimodular materials. The 
flexure problem is chosen to investigate this response since the two 
regions, one of tension and one of compression, can be identified eas- 
ily using simple intuition. Two distinct problems are considered: the 
first is where the stress is assumed continuous across the boundary of 
the two regions, which assumption has a sound physical basis. The 
second problem considered is more speculative: it is where disconti- 
nuities of stress are allowed. It is shown that such discontinuities are 
necessarily small for many applications, but might nonetheless provide 
an explanation for the damage incurred by repeated flexure. Some ex- 
perimental evidence of the possibility of bimodularity in elastomers is 
also presented. 
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1 Introduction 



Many materials seem to behave differently in tension and compression. Such 
materials are sometimes termed himodular materials. Such bimodularity has 
been observed, for example, in rocks (Lyakhovsky et al. [1]), nacre (Betoldi 
et al. |2]), and soft biological tissue, especially cartilage (Soltz and Ateshian 
[3]). Although the concepts of tension and compression seem intuitively ob- 
vious, a precise definition of these terms is a non-trivial problem that is still 
open. Perhaps the most definitive formulation of these concepts is given in 
the study by Curnier et al. [Ij, who formulate their theory in the context of 
unconstrained non-linear elasticity. The linearization of their theory yields an 
unambiguous and precise definition of tension and compression for infinites- 
imal strains, but the situation for non-linear deformations is less clear-cut, 
with the regions of tension and compression corresponding to positive and 
negative values respectively of an undefined functional of the Green-Lagrange 
strain tensor. 

This difficulty is avoided here by considering the problem of flexure, where 
intuition suggests the appropriate regions of tension and compression. If a 
rectangular bar is bent by terminal couples, with the faces along the length 
of the bar assumed to be traction-free, then the region of tension corresponds 
to the region where imaginary fibres originally aligned along the length are 
extended and the region of compression is where the fibres are contracted, 
as illustrated in Figure [T] below. This problem is, perhaps, the benchmark 
problem that all theories of bimodularity must describe; specifically, a general 
theory of when a material is in tension and compression should coincide with 
our intuitive notion of tension and compression in flexure. 

The problem of flexure is considered here within the context of non-linear, 
isotropic, incompressible elasticity, primarily because the mathematical anal- 
ysis is considerably simpler than that for unconstrained materials. First for- 
mulated and solved by Rivlin [5] , it has since been studied extensively in the 
literature (see, for example. Green and Zerna [6] and Ogden [^). There is 
continuing theoretical interest in this problem, as can be seen, for example, 
in the recent studies of Kanner and Horgan [B] and Destrade et al. fSj. Rivlin 
[S] showed that if a circular, annular sector is assumed for the deformed con- 
figuration, then an elegant solution to a natural boundary value problem 



can be found. We recall this derivation in Sections |3.1| and |3.2[ There, the 
solution is simplified by assuming that the turning angle a through which 
the beam is bent is specified, as opposed to the more usual assumption of 
specifying the moment. 

Next (^3.3), we establish the general expression for the stress difference 
across the neutral axis, which separates the region of "compression" from 
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Figure 1: In-plane section of a bent block of an isotropic, incompressible, 
non-linearly elastic solid. The bending angle is a. The inner and outer radii 
of the bent faces are and r?,, respectively. Along the arc line at r = r„, the 
circumferential line elements preserve their length during bending. In the 
region of "compression", < r < r„, they are contracted. In the region of 
"tension", r„ < r < r^, they are extended. 

the region of "tension". The consequences of assuming that the stress is 
continuous across the neutral axis are first explored. It is shown in Section 
|4] that for the case of a Mooney-Rivlin solid, or equivalently, of a general 
incompressible solid in the third-order approximation of non-linear elasticity, 
knowledge of /i, the ratio of the shear moduli in tension and compression, is 
all that is required to solve the boundary value problem in its entirety. This 
solution is then used to illustrate some possible consequences of bimodularity 
in flexure. The possibility of a stress jump across the neutral axis is also 
briefly considered with an estimate given as to its possible magnitude 

The problem considered here is one of the simplest non-linear problems in 
solids that involve structural changes induced by applied mechanical forces. 
It is hoped that this analysis can give insight into more complicated problems 
involving such changes and also can provide a reference solution for non- 
linear flexure problems where structural changes other than simple changes 
in material constants are considered. The problem of flexure is studied as 
it is very often the dominant mode of deformation, in the sense that the 
applied forces encountered are small compared to those necessary to effect 
other deformations, for many materials that experience structural changes, as 
is shown, for example, in the work of Lua et al. pi)]. It is also of importance 
in the theoretical investigation of the consequences of proposed models of 
complex behaviour in solids (see, for example, Rajagopal et al. [TT]). 
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Although the idea of elastomers behaving differently in tension and com- 
pression has an intuitive appeal, experimental evidence of this is scant, pre- 
sumably because the effect of bimodularity has not previously been consid- 
ered important. Some indirect evidence of this effect in rubbers is presented 
next. 

2 Experimental data on bimodularity 

Rubbers are traditionally modeled as being non-linearly elastic, incompress- 
ible materials. The recent data of Bechir et al. |!l2] suggest that, for at 
least some natural rubbers, there seems to a different response in tension 
and compression. 

By examining closely the data displayed on their Figure 8(a), we can es- 
timate the slope of the stress-strain curve in the linear region of compression. 
We obtain an excellent (visual) fit for the first recorded five data points by 
taking a straight line passing through the origin with slope approximatively 
equal to 8.33 MPa. For the region of tension, we use their experimental 
data (kindly provided by H. Bechir), see Table 1. Here, depending on how 
many points we estimate to be in the linear regime, we find that the slope is 
somewhere between 3.94 and 5.84 MPa. This suggests that for the natural 
rubber NR70 of Bechir et al. [12], the ratio of to E~ , the infinitesimal 
Young moduli in tension and in compression, is within the range 

0.47 < E-^/E- < 0.70. (2.1) 

Table 1. Experimental data of Bechir et al. fT^ in the early (linear) stages 
of the uniaxial tensile region: Cauchy stress (MPa) Vs elongation (m/s). 



e 


0.0514 


0.1016 


0.1518 


0.2023 


a 


0.3000 


0.4899 


0.6525 


0.7971 



We conducted similar experiments on soft translucent silicone [Feguramed 
GmbH] . Figure [2] shows representative results obtained from uni-axial com- 
pression and tension tests, in the neighborhood of the unstressed state. We 
carried out several tests in each regime, on specimens with various sizes, and 
found that in the linear tensile region, the slope was approximatively 0.45 
MPa (notice that the silicone is about 20 times softer in tension than the 
NR70 of Bechir et al.), while in the compressive region it was approxima- 
tively 0.25 MPa. This suggest that in contrast to rubber, silicone is stiffer in 
tension than in compression, with 

E+/E- ~ 1.8. (2.2) 
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Figure 2: Uni-axial compression and tension tests of translucent silicone. The 
original dimensions of the block in compression were 50 mm (depth) x 70 mm 
(width) X 50 mm (height); there, the elongation was measured as (recorded 
displacement) /(original height); the Cauchy stress was computed from the 
recorded nominal force applied by the Tinius Olsen machine; the contact 
areas between the compression platens and the specimen were generously 
lubricated. The dogbone specimen used in tension had an original section 
area of 2 mm x 6 mm; there elongation was tracked directly by LASER 
monitoring. 



Of course, we must acknowledge that external factors, other than bi- 
modularity, can explain that / ^ 1, most importantly, those due to 
experimental error and protocol. Indeed uniaxial compressive tests such as 
those presented above are harder to implement than tensile tests, because a 
perfect lubrication between the platens and the sample is required in order 
to ensure homogeneous deformations. Nonetheless we note that there exists 
standardized protocols for compressive tests of rubbers (see ISO 7743 [13], 
ASTM D575 [HI, and also Gent [15]). It seems thus perfectly legitimate to 
investigate what would be the consequences of bimodularity for elastomers. 

Finally we remark that Soltz and Ateshian [3J used bimodularity to model 
the tension-compression nonlinearity of articular cartilage. There, the mean 
values obtained from 9 specimens give E^ / E^ = 12.75/0.6 = 21.25. Simi- 
larly, Bertoldi et al. [2] recently showed that nacre (mother-of-pearl) must be 
considered orthotropic and bimodular in order to interpret correctly the avail- 
able experimental data. They estimate that for nacre, 0.5 < E^ /E~ < 0.8. 
They also carry out calculations for the plane strain bending of a orthotropic, 
bimodular beam, within the theory of linear elasticity. In what follows here, 
we focus on the large bending of isotropic, incompressible, bimodular beams 
in non-linear elasticity. 
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3 Bending 



3.1 Bending deformation 

The fundamental assumption introduced by Rivlin [S] to model the non-linear 
flexure of an incompressible beam is that a beam of length L and thickness 2 A 
is deformed under applied terminal moments into a circular, annular sector. 
For definiteness, assume that the faces X = ±A are deformed into the inner 
and outer faces of radii and r^, respectively, of the annular sector and 
that the faces Y = ±L/2 are deformed into the faces 6 = ±a, where a is a 
specified constant. The turning angle a, is restricted to lie in the range 

0<a<7r, (3.1) 

which only allows a beam to be bent into at most a circular annulus. 

Adopting a semi-inverse approach, Rivlin j5] shows that the following non- 
homogeneous deformation field is a solution to the equilibrium equations of 
non-linear incompressible elasticity: 

r = ^/2{L/a)X + D, 9 = aY/L, z = Z, (3.2) 

where {X, Y, Z) and (r, 6*, z) denote the Cartesian and cylindrical polar coor- 
dinates of a typical particle before and after deformation, respectively, and 
D is a constant. Hence, the inner and outer radii of the deformed curved 
surfaces are determined by 



r,,b = ^DT'2{L/a)A. (3.3) 

Adding and subtracting these equations then yields 

D = {rl + rl)/2, rl-rl = AAL/a. (3.4) 

The corresponding deformation gradient tensor, is given by 

-F = diag (Ai, A2, A3) = diag (a, A~\ l) , where X = L/{ar), (3.5) 

denoting the principal stretches by Ai, A2, A3. Here A2 = A^^ is the stretch 
experienced by circumferential line elements: when A2 > 1 (A < 1), the 
line elements are extended by the bending; when A2 < 1 (A > 1), they 
are contracted. Also, A3 = 1 at all times, showing that Rivlin's solution 



(3.2) is a plane strain deformation. Note that A1A2A3 = 1, showing that the 
deformation respects the internal constraint of incompressibility. Finally, it 
is clear that F corresponds to a non-homogeneous deformation because the 
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radial stretch Ai varies as the inverse of r. In particular, it takes the following 
values 

Xa = L/ (ava), \b = L/ (arfc), (3.6) 
on the inner and outer faces of the bent beam, respectively. These notations 



allow us to rewrite (3.4)2 in non-dimensional form as 

X,' - Xf = 26, (3.7) 
where e is the product of the beam aspect ratio by the turning angle [9]: 

e = {2A/L)a. (3.8) 

3.2 Homogeneous (unimodular) beams 

For homogeneous, incompressible, elastic materials, the corresponding prin- 
cipal Cauchy stresses are 

Trr = -p + XiW,i, Tee = -p + X2W,2, (3.9) 

where p is an arbitrary scalar field, W = W{Xi, X2, X^) is the strain-energy 
function and the comma subscript denotes partial differentiation with respect 
to the appropriate principal stretch. The equations of equilibrium determine 

p as 

p = j (XiW^i - X2W,2) r~Mr + Ail^,! + K, (3.10) 
where K is an arbitrary constant. It therefore follows immediately that 

Trr = J (XiW^i - X2W,2) r^Mr + K, Tee = Trr + AsVT.s - Ail^.i. (3.11) 

Now define the function W{X) as 

W{X) = W{X,X~\1), (3.12) 

which is assumed henceforth to be a convex function. Then 

XW' = Aiiy,i - X2W^2, (3.13) 

where the prime denotes differentiation. The stress components can then be 
written simply as functions of A as 

T„ = W + K, Tee = Trr - XW' . (3.14) 
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The curved surfaces of the bent beam are assumed to be free of traction. 
This assumption then yields 



K 



-W{Xa) 



W{Xb) = W{Xa). 



(3.15) 



Only isotropic materials are considered in this paper. For these materials, 
W{Xi, A2, 1) = W{X2, Ai, 1) , or equivalently, 



W{X) = w{x-^). 



Then (3.15)2 yields 



W{X,) = W{X^) = W{X,') = w{x-') 



(3.16) 



(3.17) 



There are two obvious solutions to these equations: the first is Xa = X^, 
which by (3.6) gives = Vf,, a physically unacceptable solution; the second 
is Xa 



X 



, which gives 

XaXb 



or 



2 r 2 

a ran = L . 



(3.18) 



In that latter case, (3.6) yields a quadratic equation for A^, with the following 
unique physically acceptable solution: 

Xa = 



or 




4v42 



-2A 



(3.19) 



which completely determines the deformed configuration. Unusually, it is 
independent of the form of the strain- energy function (provided that (3.18), 
which is sufficient for (3.17) to be satisfied, is also necessary). From the 
specialization to isotropic materials follows also immediately from (3.13) that 



0. 



(3.20) 



Hereafter we consider the boundary value problem where equal and op- 
posite moments are applied to the ends of the beam at F = ±L/2, whilst 
the inner and outer curved faces of the beam are traction-free. 



3.3 Bimodular beams 

The neutral axis in flexure is the circumferential material line (originally 
parallel to the major axis of the beam in the reference conflguration) with 
unchanged length, see Figure [1} The neutral axis is therefore given by the 
radius r„ in the deformed conflguration such that 

Tn = -, (3.21) 
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assuming that 



ra<rn< n- 



(3.22) 



The neutral axis is the natural, intuitive, boundary in flexure between the 
regions in "tension" and "compression", defined by r > r„ and r < r„, 
respectively. 

Assume now that a material behaves differently in these two regions. 
The turning angle a is assumed to be constant across these two regions. 
Assuming also that the radial deformation field is continuous across r„ then 
means that both the deformation field and the deformation gradient tensor 
have the same form in both regions. Denote quantities associated with the 
tensile and compressive regions by the superscripts '+' and respectively. 
It follows from (3.14) that the radial stress therefore has the form 



T+ = W+{X) + K+, T- = W-{\) + K- 



(3.23) 



where K^,K~ are constants and A is defined by (3.5). Satisfying the stress 
free boundary conditions on the upper and lower curved surfaces, and solving 
for K~^, K~ then yields 



T+ = l^+(A)-W^+(Ab), 
Along the neutral axis then 



W-{\)-W-{\a). 



(3.24) 



W+{1) - W+{\), T-{r^) = W-{1) - W-{Xa). (3.25) 



The usual assumption of zero strain energy in the reference configuration 
yields ^ ^ 

iy+(l) = ly-(l) = 0, (3.26) 

and so we obtain the following relation for the normal stress difference across 
the neutral axis: 



AT,. 



(3.27) 



It follows from (3.11) and (3.20) that the difference in the hoop stress is 

- T^,{r^) = W-{K) - W^{\). (3.28) 



exactly the same: 
ATf 



T+( 



We call AT this stress difference: AT = AT,., = ATeg. 

Consider now that there is no stress discontinuity across the neutral ra- 
dius. Then AT = and 



W-{Xa) = W+{Xb). 



(3.29) 
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Recalling (3.7), this equation can be rewritten as the defining equation for 



Xa as follows: 



w~iXa) = w+iXaii + 2eXly 



(3.30) 



With Xa, necessarily > 1, determined in this way, the deformed configuration 



is completely determined, 
problem solved by Rivlin (see (3.19)) 



Note that, in contrast to the standard fiexure 
the deformed configuration depends 



here on the form of the strain-energy function. One such form is considered 
in the next section. 



4 An example 

In this section we solve the fiexure problem in the case where the solid is 
modeled by the Mooney-Rivlin strain energy density: 

W='^{l + f){Xl + Xl + Xl-3) + |(1 - mlXl + A^A^ + XlXl - 3), (4.1) 

where /i and / are constants, or equivalently [I7j, by the third-order expan- 
sion of W in the weakly non-linear elasticity approximation flEl [T9] : 

W = i2ti{E^) + '^tT{E^), (4.2) 

where fi is the shear of modulus of second-order elasticity, ^ is a third-order 
non-linear Landau constant, and E is the Green-Lagrange strain tensor (with 
eigenvalues (A^ — l)/2). 

In both cases we find that 

WiX) = I (A' + A-2 - 2) . (4.3) 



Then (3.30) reduces to the following cubic: 

2ex^ + [l - 4e - /i (1 - 2ef] a;^ - 2 [1 - e - /i (1 - 2e)] s + 1 - /i = 0, (4.4) 



where x = A„, e is defined in (3.8), and p, = n'^ / is the ratio of the 
shear modulus in the region of tension by the shear modulus in the region of 
compression. 

Focusing now on bimodular solids which are stiffer in compression than 
in tension, we take < /t < 1. Calling f{x) the cubic on the left hand 
side of ( |0| ), we find that /(I) = -4/ie^ < and /(oo) = 1 - /i > 0, 
ensuring that there always exists a relevant root for Xa- Setting fi = 1 
recovers the case where there is no difference in material properties in tension 
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and compression, and (4.4) then yields (3.19). Letting /i — )■ in (4.4), which 



physically corresponds to the case where /i+ << yields = 1: in other 
words, an infinitely stiff bar cannot be bent, as expected intuitively. Figure 
|4] displays several graphs of Xa as a function of e, for different values of the 
parameter fi. It is clear that the effect of the tension-compression difference 




X8 



Figure 3: Variations of A^, the radial stretch ratio on the inner curved face, 
with e, the product of the bar aspect ratio by the turning angle, in the case of 
a bimodular Mooney-Rivlin solid. The ratio of the tension shear modulus to 
the compression shear modulus takes the values: /t = 0.0 (horizontal dotted 
line), 1/20, 1/10, 1/4, 1/2, 1/1.1 (full lines), and 1.0 (other dotted line). 

becomes more pronounced with increasing angle and/or aspect ratio. 



5 Possibility of a stress discontinuity 

Finally, the possibility of a stress discontinuity across the neutral axis is con- 
sidered. The size of the stress discontinuity cannot be determined without a 
further assumption. It will be assumed here that the deformed configuration 
for the bimodular material is the same as that in the classical problem of 
fiexure, 



As shown in Section 3.2, the deformation of the classical fiexure prob- 
lem is independent of the form of the strain-energy function and therefore, 
in particular, independent of which of the tension/compression forms is to 
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assumed. For the classical problem, 

Xt = X-\ (5.1) 



according to (3.4) and, from (3.27), the stress discontinuity in the bimodular 
beam therefore has the value 

AT = W-{\a)-W+i\-'). (5.2) 



Noting (3.16), this can be simplified as 

AT = W-{Xa)-W+{Xa). (5.3) 

Now the material is isotropic, and without loss of generality, the strain- 
energy function can be assumed to have the form 

W{X) = W{I), where J = + A'^ - 2, (5.4) 

because in bending, the first and second principal invariants of the Cauchy- 
Green strain tensor are both equal to /. The problem of flexure is a typical 
example of a deformation where the strain, quantified here by /, is small but 
with possibly mo^rate or large deformations. Therefore the general strain- 
energy function W{X) can be closely approximated by the linearisation of 



W{I), i.e., by the form (4.3) of the strain-energy function. It follows that 
AT can be closely approximated in general by 

AT = - (A^ + - 2) = 2 (/i- - (v/^^ - l) . (5.5) 

Hence, in general, it can be seen that the stress discontinuity cannot be 
large because the difference in the shear moduli is multiplied by a small 
term. We would like to propose this stress difference as a possible cause of 
damage in the flexure of hyperelastic materials. We note that simply flexing a 
hyperelastic bar is not likely to cause much damage but, over time, flexing and 
unflexing the same bar is likely to cause the accumulation of small amounts 
of damage due to the stress discontinuity, with inevitable consequences for 
the integrity of the component. 
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